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ABSTRACT 

The common-envelope (CE) phase is an important stage in the evolution of 
binary stellar populations. The most common way to compute the change in 
orbital period during a CE is to relate the binding energy of the envelope of the 
Roche-lobe filling giant to the change in orbital energy. Especially in population- 
synthesis codes, where the evolution of millions of stars must be computed and 
detailed evolutionary models are too expensive computationally, simple approxi- 
mations are made for the envelope binding energy. In this study, we present accu- 
rate analytic prescriptions based on detailed stellar-evolution models that provide 
the envelope binding energy for giants with metallicities between Z = 10 -4 and 
Z = 0.03 and masses between 0.8 M Q and 100 M , as a function of the metallicity, 
mass, radius and evolutionary phase of the star. Our results are also presented in 
the form of electronic data tables and Fortran routines that use them. We find 
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that the accuracy of our fits is better than 15% for 90% of our model data points 
in all cases, and better than 10% for 90% of our data points in all cases except the 
asymptotic giant branches for three of the six metallicities we consider. For very 
massive stars (M > 50 M ), when stars lose more than ~ 20% of their initial 
mass due to stellar winds, our fits do not describe the models as accurately Our 
results are more widely applicable — covering wider ranges of metallicity and 
mass — and are of higher accuracy than those of previous studies. 

Subject headings: stars: evolution, stars: fundamental parameters, stars: mass 
loss, binaries: close 



1. Introduction 



The common-envelope (CE) phase ( jPaczynskilll976uWebbinklll984l ; iTaam fc Sandquist 



20001 ) is an important event in the evolution of many binaries and is used to explain large 
numbers of observed compact binaries, such as X-ray binarie s, cataclysmic variables (CVs) 
and double degenerates, as well as bipolar planetary nebulae ( Ide Koollll990[ ). The CE phase 
is believed to be initiated when an evolved, giant star is in orbit with a more compact 
companion and fills its Roche lobe. When such a donor star has a deep convective enve- 
l ope and/or the mas s ratio of the system is sufficiently large, unstable mass transfer occurs 
( IRasio fc Liviolll996l ). resulting in a fast expanding envelope which quickly engulfs the com- 
panion. Inside this common envelope, the binary orbit of the core of the giant and the 
secondary star shrinks due to friction and torques, leading to either a compact binary con- 
sisting of these two objects or their merger. The energy generated through orbital shrinkage 
heats and eventually expels the envelope. The whole process is thought to occur on time 
scales much shorter than the nuclear-evolution timescales of stars (< 10 3 yr) so that the mass 
of th e giant's core does not change, and the secondary is not affected (ITaam fc Sandquist 

2oooh . 



Despite the fact that CEs are an important ingredient for modelling stellar systems 
and populations, the details of the process are poorly understood. Three-dimensional hy- 
drodynamical models have provided detailed simulations for the first month or so of CE 
evolution, in which significant orbital shrinkage takes place, but due to the large differences 
in scale between a red-giant envelope and its core, the precise outcome of the CE phase 
cannot yet be predicted (e.g., ITaam &: Rickerll2006l ). Instead, a very simplified scenario is 
commonly used, where the binding energy is equated to the difference in orbital ener gy be- 
fore and after the CE in order to predict the post-CE orbital period ( jWebbinkl Il984l ) . Part 
of the uncertainty in CE evolution is allowed for in the CE parameter acE, which speci- 



-3 - 



fies the efficiency with which orbital energy is used to expel the envelope. Especially for 
population-synthesis codes, where the evolution of millions of stars must be computed and 
models are very basic, the stellar structure needed to compute the envelope binding energy 
is not available and the envelope-structure parameter A cnv (see Eq.E] and the discussion in 
Sect. 15.4) ) is used to approximate the binding energy from basic stellar param eters. In many 



0.5 



e.a 



de Kool et al. 



studies, this parameter has been assumed const ant, typically A cn 

1987 ; iNelemans et al.ll2000l ; iHurley et al.l 120021 ). or acE^env = 1-0 or 0.5 (iBelczynski et al. 



20081 ) , whereas in reality it can vary wi ldly during the evolu t ion of a star, especially on the 



asy mptotic giant branch 



(AGB) {e.g., iDewi fc Taurisl l2000l ; Ivan der Sluvs et all 120061 ). In 



fact, Ivan der Sluys et al.l (120101 ) show that the latter assumption implies «ce > 1 for about 
60% of the CEs that occur in a stellar population of solar metallicity. No te that an alter- 



nativ e scenario for envelope ejection, based on angular-m omentum balance (INelemans et al. 



2000 ; Nelemans fc Tout 2005 ; van der Sluvs et~al] bood ). does not depend on the envelope 



binding energy and is therefore not affected by our results. 



Yakut fc Eggletonll2005l . and 



In this paper, we use the stellar-evolution code ev (e.g., 
references therein) to compute detailed stellar-evolution models for a range of masses (0.8 — 
100 M ) and six different metallicities, and compute their envelope binding energies through- 
out their evolution. We fit the binding energies as a function of basic stellar parameters 
(metallicity, mass and radius) and provide the results of these fits, as well as routines to 
compute the envelope binding energy for any combination of mass and radius, for a given 
metallicity in our grid. We find that our analytic prescriptions describe 90% of our model 
data poi n ts wit h an accuracy better than 10-15% for the metallicities provided. Recently, 
Xu fc Lil (120101 ) have used the same approach, providing fits for 14 discrete masses and two 
different metallicities. Our study improves on that by using the stellar mass as a fitting 
parameter, so that any mass — and masses from a wider range — can be used, and by 
providing a larger number and wider range of metallicities. In addition, we show that our 
fits are largely independent of the wind mass loss, allowing for many wind prescriptions, and 
we provide a more detailed analysis of our accuracies. 

In Section [2] of this paper, we present the stellar-evolution code used and the models 
generated with it. Section [3] describes the method we use to fit the binding energies and in 
Section H] we present our results and their accuracies. Section \5\ contains a discussion of this 
study and we present a summary and conclusions in Section [6j 
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2. Stellar evolution 

Here we describe the stellar-evolution code that we use for this study, the relevant details 
of the physics involved, as well as the actual grids of detailed stellar-evolution models that 
are computed to perform our fits and determine their accuracy. 



2.1. Stellar-evolution code 



We compute our stellar models using a version of the binary ste 



(also known as STARS or TWINjj, developed by Eggleton (lEggleton 



20051 . and references therein) and updated as described in iPols et al. 



1971 



lar evolution code ev 



1972 



Yakut fc Eggleton 



(119951 1 . In the code, 



convective mixing is modelled by a diffusion equation with a ratio of mixing length to pres- 
sure scale height of l/H p = 2.0. Convective overshooting in the core is taken into ac- 
count on the main sequence (MS) for stars with M > 1.2 M Q and on the horizontal branch 
(HB) for stars of all masses. On the MS, we use an overshooting parameter 5 ov = 0.12 for 
stars with M > 2.0 M , which corresponds to an overshooting length of about 0.3 H p , and 
0.0 < <5 ov < 0.12 increasing linearly with mass for stars with 1.2 M < M < 2.0 M Q . On the 
HB, we use 5 ov = 0.12 for all models. The code cannot evolve a model through the helium 
flash, the violent ignition of helium in a low-mass star with a degenerate core, but it auto- 
matically replaces the model at the moment of helium ignition with a tailored model of the 
same total mass and core mass in which helium has just ignited. For stars above a certain 
mass (M > 2.1 M Q for Z = 0.02), helium ignites non-degenerately and this intervention is 
not needed. 

We define the helium-core mass M c as the mass coordinate below which the hydrogen 
abundance does not exceed 10%. We compute the binding energy of the hydrogen-rich enve- 
lope, -Ebind; by integrating the gravitational and internal energies over the mass coordinate 
of the model, from the core-envelope boundary to the surface of the star M s : 



bind 



Mb Gm 

M c 



T 



r\m) 



The term E- mt is the internal energy per unit of mass, which contains terms such as the 
thermal energy of the gas and the radiation energy, but not the recombination energy. More 
details regarding these assumptions are provided in Ivan der Sluys et al.l (120061 ). 



lr The current version of ev is obtainable on request from eggletonl@llnl.gov, along with data files and 
a user manual. 
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In this study, we compute grids of models with six different metallicities: Z = 10 4 , 
Z = 0.001, Z = 0.010, Z = 0.015, Z = 0.020, and Z = 0.030. The initial hydrogen 
and helium abundances (X and Y, respectively) of our model stars are a function of the 
metallicity Z, and are given by X = 0.76 - 3.0Z and Y = 0.24 + 2.0Z. 

For three of these metallicities, Z = 10~ 4 , Z = 0.020, and Z = 0.030, we calculate 
additional grids tha t inclu de mass loss via stellar winds. We use a prescription that was 



inspired by iReimersl ( 119751 ) : 



3.16- lO^Moyr" 1 (m) (l) (jfa<^. 
M = - V -mm{ ° W W W , (2) 

g.Gl-lO-^Moyr- 1 (±) 

where we set t] = 0.2 for this study. This wind prescription dominates the stellar winds for 
lower-mass stars on the giant branches, and the upper prescription in Eq.|2]is used in most of 
those models, except where l-EWndl is small, e.g., for stars near the tip of the as ymptotic giant 



branch (AGB). For massive, luminous stars, we use the wind prescription by Ide Jager et al. 



(11988( 1 . 



2.2. Stellar-evolution models 

We computed several grids of single-star models for a range of metallicities and for each 
metallicity, for a range of initial masses. For these grids, we assume that there is no mass 
loss due to stellar winds. However, for three metallicities {Z = 0.02 ("solar") and the two 
extremes Z = 10~ 4 and Z = 0.03), we computed a grid of models where wind mass loss is 
included, in order to determine how the accuracy of our fits is affected by stellar winds and 
how we can correct for that (see Sect. 13.41 and 14. 3|) . 

We selected 73 zero-age main-sequence (ZAMS) masses between 0.8 M & and 100 M© in 
each grid, distributed uniformly in log M: 

log Mi ~ -0.038660 + 0.029124 • i (3) 

resulting in M ~ 0.80, 0.86, 0.91, 93.5, 100.0 M Q . We ignore stars that have a main- 
sequence lifetime longer than 13 Gyr for a given metallicity. We use the models with an 
odd value of % to perform the fits, whereas the models with an even value for i are used as 
verification, to determine the accuracies of our fits for models that they are independent of. 

The full grids of detailed stellar-evolution models are computed for all metallicities 
considered. For the three stellar-wind grids, only odd-i models are calculated, since this is 
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all that is required to determine the mass- loss correction factor described in Sect. l3.4l A 
Hertzsprung-Russel diagram (HRD) for selected stellar models from our grid with Z = 0.02 
is shown in Fig.[TJ 



3. Method of fitting 

With the data generated as explained above we construct fit functions to model the 
behavior of the binding energy. In this section, we describe the procedure by which these 
fits are obtained. Section 13.11 describes our choice of fitting parameters, in Sect. 13. 21 we 
describe how we divide our evolutionary tracks into groups which will be fitted separately, 
and Sect. 13.31 contains the formulation of the fitting functions and a description of the fitting 
procedure itself. Finally, in Sect. 13.41 we present a correction factor that takes into account 
the most important effects of wind mass loss. 



3.1. Choice of fitting parameters 

After some experimentation we find that the binding energy is best described by the 
stellar mass and radius, as opposed to other basic parameters like age or core mass. Such 
a relationship of the binding energy with the mass and radius of a star was also found by 



Dewi fc Taurisl ( 120001 ) . Because the ranges of E hind , R and M may span several orders of 



magnitude, we use the (base-10) logarithm of these variables. 

We do our fits separately for each metallicity and do not provide fits that take into 
account metallicity as a parameter. The reason for this choice is that with each added fitting 
parameter, subtle features in the models that were well-described by the other parameters 
are washed out when the new dimension is added. However, we compare the changes due 
to the metallicity in the log -EWnd — log R plane for models with a given mass, and find that 
these changes are relatively smooth, see Fig.[3j This suggests that the user could compute 
the envelope binding energy for the two metallicities that bracket the desired value of Z and 
interpolate linearly in log -Ebmd i n order to estimate the binding energy for the given star. In 
practice, however, this grid may be sufficiently dense that the user can take the value of Z 
that is closest to the desired value. 
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Fig. 1. — Hertzsprung-Russel diagram for selected models with masses of 0.91, 1.36, 2.04 
3.05, 4.57, 6.84, 10.22, 15.3, 22.9 and 34.2 M , from our Z = 0.02 grid. The different line 
styles and colours indicate different evolutionary stages. Dotted (grey) lines are used for 
the main sequence, (part of the) Hertzsprung gap and horizontal branch, phases where the 
star either has no deep convective envelope or cannot fill its Roche lobe. The solid (red) 
lines indicate the red-giant branch (LMR), the dashed (blue) lines show the asymptotic giant 
branch (LMA) and the dash-dotted (magenta) lines indicate the high-mass models (HM), 
which have only one giant branch. 
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3.2. Dividing the evolutionary tracks 

A CE phase can only be initiated by a star that is at its maximum radius so far in 
its evolution. Because of this we consider only evolutionary phases in which a star expands 
(almost) monotonously. This results in the division of the evolutionary tracks into four 
different groups, and we develop fits for each group separately. 

First, we divide our set of stellar-evolution models into low- mass and a high- mass groups. 
The dividing mass between the two groups lies around 12 M , depends lightly on metallicity, 
and is given in Table [TJ 

The low-mass stars in our grid (M < 12 M ) evolve through a clear horizontal-branch 
(HB) phase between the red-giant branch (RGB) and asymptotic giant branch (AGB). In 
this phase, the star burns helium in its core and the star's radius is smaller than it was on the 
RGB, so that this phase can be ignored for common-envelope evolution. In addition, most 
of these low-mass stars undergo a first dredge-up on the RGB, during which the hydrogen- 
burning shell meets less-processed material and the star as a whole shrinks slightly. We define 
the beginning of the RGB as the moment where the hydrogen abundance in the stellar core 
becomes zero. As a consequence, the envelope binding energy in the Hertzsprung gap (HG) 
is also described by the RGB fits for most stars (see Fig. CP . We define the end of the RGB 
as the point of maximum radius before the core helium abundance drops below half of its 
absolute maximum value over the stellar lifetime. The beginning of the AGB is taken to be 
the point where the core helium abundance drops below approximately 10% (we vary this 
value slightly for different stars). We find that this roughly defines the point where the stellar 
radius begins to increase again after the HB. For all of our models, the "AGB" begins well 
before the radius surpasses its maximum on the RGB, so that we cover more than the range 
of radii for which a CE can be initiated. As the tip of the AGB, we take the model where 
the star reaches its maximum radius in its entire evolution. We find that the accuracies of 

Table 1: Dividing mass between the low-mass (LMR/LMA) and high-mass (HM) models 
for each metallicity considered 





Z = 


10" 4 


11.7 M 


z = 


0.001 


11.7 M 


z = 


0.01 


10.2 M 


z = 


0.015 


11.7 M Q 


z = 


0.02 


11.7 M q 


z = 


0.03 


13.4 M 
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our fits are greatly improved if we treat the tracks before and after this dredge-up phase 
separately. Since the dredge-up occurs at different radii for models with different masses, we 
include a simple function to describe this radius as a function of the stellar mass: 

4 

-Rlmr = o»; 

with the coefficients Oj given in Table [2j 

As a result, we split the tracks from the low-mass models into three groups, LMR1 
for the low-mass early RGB (before the dredge-up), LMR2 for the low-mass late RGB and 
LMA for the low- mass AGB. 

For the models with higher masses (M > 12 M & ) there is no HB and no first dredge-up, 
and we keep these tracks intact and put them in a fourth group, which we call HM. Hence, 
these tracks are defined between the moment hydrogen is exhausted in the core near the 
TAMS and the maximum radius ever achieved in the evolution of each star. The lowest- 
mass model (whose mass is ~ 12M Q , see Table [I]) included in this high- mass group is also 
included in the low-mass RGB and AGB groups to ensure continuity between the regions (we 
use the phenomenological parallels in the log-Ebind — logi? curves (see Fig. [2]) of lower-mass 
models to determine where to divide this track between RGB and AGB). 

We interpolate the evolutionary tracks from each stellar model for each group in the 
log-Ebind — log-R plane, in order to give each model the same number (200) of data points, 
which gives them an equal weight when fitting along the mass axis, and with a constant 
density along their tracks in this plane, which prevents biasing the fits towards particular 
regions in these curves. The only exception are the groups LMR1 and LMR2, which are 
given 200 data points together. Figure [2] shows the model tracks in the LMR, LMA and HM 
groups for a selection of the models in our Z = 0.02 grid. 



MS 



Table 2: Coefficients for the fit of the division between the groups LMR1 and LMR2 in 
Eq.SJ 
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a i 


a 2 


«3 


a4 


10" 4 


0.604723 


1.39724 


-0.0824963 


1.14143 





0.001 


0.4566 


1.1866 


2.39088 


-3.05404 


1.4049 


0.01 


0.282174 


1.14938 


1.88445 


-1.0823 





0.015 


0.251818 


1.21049 


1.6349 


-0.83691 





0.02 


0.240637 


1.08922 


1.95318 


-1.03218 





0.03 


0.234888 


0.897294 


2.51995 


-1.41411 
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3.3. Selecting the fitting functions 



As described above, we choose to find a description for logi?bind as a function of logM 
and logi?. As a general functional form, we select the polynomial of an as yet undefined 
order in both the instantaneous mass M and radius R. Hence, the general form of our fitting 
function is 



E 



bind 



log , 
A(M ,M)-J2 



E + 



a, 



log 



M 



log 



R 

Rq 



(5) 



with Eq = 33.29866, a m , r the fitting coefficients and m and r integer indices and exponents. 
The factor A depends on the initial (ZAMS) and instantaneous mass and corrects for effects 
due to wind mass loss in high-mass stars. It is of course equal to unity for stars that do not 
experience (significant) mass loss and is discussed in more detail in Sect. l3.4l 

Finding a good and accurate prescription for the binding energy then consists of two 
parts. First, we must determine the orders of each of the polynomials in Eq.[5j In fact, 
we choose to use a liberal definition of the word polynomial, and rather than selecting the 
orders, we vary the ranges of m and r for which a m r has non-zero values in order to find the 
narrowest ranges that still describe i?bind with good accuracy, possibly including negative 
values for m and r. Second, for each of our choices of the ranges of m and r, we must find 
the optimal set of fitting coefficients a m ^ r and quantify the accuracy of the fit, in order to 
compare to fits that use different polynomials. 

For the actual fitting of the data points, we use only the odd-numbered mass models 
from the grid described in Sect. 12.21 (the masses with odd values of i in Eq.[3]). For each 
selection of the range of m and r, we fit the data points, i.e. determine the best set of a m)T , 
using the Y 2 -minimization-ba sed fitting procedures available in Wolfram Mathematica 7 
( jWolfram Research. Inc.ll2008l ). Instead of using Mathematica's "goodness-of- fit parameter" , 
we choose three criteria to compare the accuracies of the different fits in a quantitative way. 

First, we use the newly generated fit function to predict the binding energies for all mass 
models and all data points that are used for the fit (i.e. the odd-numbered mass models), 
and determine for which fraction of these data points the following condition holds 



E 



bind, fit 



Ew, 



bind,mdl 



E 



bind,mdl 



(6) 



where we choose the value / = 0.1 for our criterion (note that here we do not use the 
logarithm of -Ebind)- The fraction of data points that fulfills this criterion for / = 0.1 is 
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called A 10 %. More specifically, we demand that for at least 90% of the data points, the fit 
has an accuracy of 10% or better; in other words, A 10 % > 90%. 

Second, we use these fits to compute the envelope binding energies for the even-numbered 
mass models in our grid, which were not used to create the fits. Again, we use the condition 
in Eq.[6]for / = 0.1 to determine the accuracy of the fit, and demand that the results are 
similar to those from the first criterion. 

Third, we select the simplest polynomial, i.e. the narrowest ranges for m and r in Eq.0 
that still pass the first two criteria. 

For some groups of data these criteria turn out to be too strict, i.e. there are no reason- 
able ranges of m and r that could meet the first criterion while still being sufficiently well 
behaved to also meet the second. For these groups we are satisfied if the fit predicts at least 
90% of the data points within 15% of the model value, i.e. we increase / from 0.1 to 0.15 for 
these groups and demand that A 15 % > 0.9. For this reason we report both A 10 % and A 15 % 
for all groups in the next section. 



So far, we have considered conservative stellar-evolution models only. However, we find 
that when stars lose a significant amount of mass due to stellar winds (i.e. stars with masses 
> 20 — 30 M Q ), our fits lose some of their accuracy. Since wind mass- loss rates are uncertain, 
and because the assumption of evolution without mass loss is unreasonable for the more 
massive stars in our grids, we introduced the correction factor A in Eq.[5J The factor is 
based on a comparison of the fit using conservative models, as discussed in the previous 
section, to grids of models which experience mass loss as described in Sect. l2.ll For models 
in the low-mass grid (LMR/LMA), this correction is not necessary, and we use A = 1. For 
the high-mass models (HM), we find that a correction factor based on the relative amount 
of mass lost since the zero-age main sequence (ZAMS) gives a reasonably good prescription 
for most models, restoring the accuracy found for conservative models (see Fig. [5]): 



where Mq is the ZAMS mass and M is the instantaneous mass of the star. Note that Ahm 
is independent of metallicity and that it is valid across the range of metallicities considered 
here. 

The fitting procedure described in this section was carried out for each of the four groups 
mentioned in Sect. l3.2l and each of the six metallicities listed in Sect. l2.2l The results of our 



3.4. Compensating for stellar mass loss 
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fits, their accuracies, and their validity for the case of moderate wind mass loss are discussed 
in Sect.SJ 



4. Results 

In this section we describe the accuracies of our fits. In section 14.11 we present the 
optimal fit parameters to describe the binding energy as a function of basic stellar parameters. 
Section 14.21 presents the accuracies of these fits and in Section 14.31 we discuss the effect of 
stellar winds on this accuracy. In Section 14.41 we discuss the domain in which our fits are 
valid. 



4.1. Analytic prescriptions 

We carried out the fitting of our data points to the fitting function in Eq.[5]as described 
in the Section [3j Here, we present the polynomials that best describe the data points, i.e., 
the ranges of m and r we used in Eq.El and the values for a m; that best describe the data 
for those polynomials. We do this for each of the four groups of our data set: LMR1, LMR2, 
LMA and HM, as described in Sect. 13.21 and defined by Tables [1] and [2] and the equations in 
the appendix and for each metallicity we consider in this study. 

The ranges for m, r are variable, since for some groups of models, the binding energy 
can be described in fewer terms than for others. The extreme values used are and 20 for 
m and —5 and 20 for r. We provide the table of coefficients in electronic form, listing m, 
r and a m>r for all cases where a m ^ r is non-zero, as well as the contents of Tables [1] and [2j 
In addition, we provide Fortran routines which can read these data files and compute the 
envelope binding energy as a function of me t allicit y, mass, radius, and evolutionary phase of 



the star (RGB or AGB) in lLoveridge et all (12010( 1 



Figure [2] shows a comparison between the detailed stellar-evolution models and fits 
for selected models with a range of masses from our Z = 0.02 grid. Figure [3] shows the 
dependence of the envelope binding energy on metallicity, by comparing models with different 
metallicities for three different masses. 
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j 1 1 1 — i — i i i i | 1 1 1 — i — i i i i | 1 1 1 — i — i i i i | r 




I I I I I ''''I I I I I ''''I I I I I ''''I i L 

1 10 100 1000 

R (R ) 

Fig. 2. — Envelope binding energy as a function of stellar radius, for a selection of models 
with masses of 0.91, 1.36, 2.04 3.05, 4.57, 6.84, 10.22, 15.3, 22.9 and 34.2 M , and Z = 0.02. 
The line styles and colours indicate the same evolutionary phases as in Fig. [Hand are used 
for the results obtained by using our fit. The black dotted lines show the original stellar- 
evolution models, which overlap with the fits in most places. RGB and AGB phases are 
disconnected here, and the lowest-mass and the three highest-mass models do not have an 
AGB phase. 
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Fig. 3. — Envelope binding energy as a function of stellar radius, for models with masses of 
M = 1.79, 13.4 and 26.2 M Q and for metallicities of Z = 1(T 4 , 10" 3 and 0.02. The line styles 
and colours indicate the value of Z as indicated in the figure. Coloured lines represent fitting 
results, while black dotted lines show the detailed models. The four groups of lines are from 
lower left to upper right: 1.79 M RGB, 1.79 M & AGB, 13.4 M and 26.2 M Q . Within each 
group, stars with a lower metallicity have a higher binding energy, except near the end of 
the AGB or high-mass tracks. 
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4.2. Accuracy of the results 

Table [3] lists the accuracies of our fits for the envelope binding energy for each of the 
metallicities considered, and for each of the model groups defined in Sect. l3.2l To express 
the accuracy of each fit, we list the percentages of the fitted data points that fall within 10% 
and 15% of the model values (A 10 % and A 15 %). 

We separately list the accuracies for the odd-numbered stellar-evolution models, which 
were used to produce the fits, and the even-numbered models, which have masses that lie 
between those of the odd-numbered models and are used only for verification. Note that the 
percentages refer to the actual (absolute) value of the binding energy, not its logarithm. The 
same results are presented in Fig.lH as cumulative histograms of the fraction of data points 
that lie outside a given accuracy. 

As these data demonstrate, all of our fits meet the criteria outlined in the previous 
section, for / = 0.15 in Eq.[HJ In fact, all but three of the fits also meet the stricter criteria 
for / = 0.1. For Z < 0.01, the features in the LMA models vary irregularly from one mass 
model to the next, so that only for our low-metallicity AGB models, the / = 0.1 criteria are 



not satisfied. iPolitano et al.l ( 120101 ) show that, for the cases where a common envelope (CE) 
leads to a merger, only 17% of the CEs are initiated on the AGB, against 83% on the RGB, 
so that the lower accuracy for our LMA results probably has a somewhat smaller impact 
than suggested by the numbers in the table alone. 



4.3. The effect of stellar mass loss 

The accuracies in the previous section apply to the prediction of the binding energies of 
stars that have not suffered mass loss due to stellar winds. Here, we compare the outcome of 
our fits, based on conservative models, to grids of stellar models that undergo wind mass loss 
as described in Sect. 12. ll Table H] lists the accuracies for this comparison for the two extreme 
metallicities in our grid (Z = 10~ 4 and 0.03) as well as "solar" metallicity (Z = 0.02), 
expressed as A 10 % and A 15 % (see Sect. l3.3p . We see that when the factor A in Eq.[5] is 
ignored, the LMR grid is hardly affected compared to Table El whereas the accuracies in the 
LMA and especially the HM grids suffer appreciably. As described in Sect. l3.4[ the correction 
factor A defined in Eq. [7] improves the accuracies for the high-mass models significantly. This 
improvement is also shown in Fig.[5j For the LMA group, we found no simple correction 
factor that is independent of metallicity. Since the drop in accuracy due to winds is not as 
large as for the LM models, we decided to leave these fits uncorrected. 

The advantage of the A factor is that our prescription for the binding energy can be 
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Table 3: Accuracies of our fits, for each metallicity considered and each group defined in 
Sect. 13.21 For an explanation of the A's, see Sect . 13. 3 1 and Eq.EJ 

Z group A io% A 15% 







odd 


even 


odd 


even 


icr 4 
irr 4 

1(T 4 


LMR 
LMA 
HM 


97.6% 
86.6% 
94.8% 


97.7% 
86.4% 
93.6% 


98.9% 
94.8% 
96.4% 


99.1% 
94.8% 
94.6% 



icr 3 


LMR 


10" 3 


LMA 


irr 3 


HM 



92.0% 90.4% 
85.3% 84.3% 
91.1% 90.0% 



95.0% 93.6% 
95.2% 92.5% 
95.5% 93.4% 



0.01 LMR 
0.01 LMA 
0.01 HM 



95.7% 93.7% 
86.3% 87.2% 
92.4% 90.0% 



97.1% 96.5% 
95.9% 95.0% 
96.2% 96.1% 



0.015 LMR 
0.015 LMA 
0.015 HM 



95.0% 95.9% 
91.4% 90.6% 
90.0% 90.3% 



96.6% 97.0% 
97.8% 96.4% 
95.5% 95.8% 



0.02 LMR 
0.02 LMA 
0.02 HM 



94.3% 92.7% 
97.1% 91.9% 
92.0% 91.7% 



96.4% 96.4% 
99.3% 96.9% 
96.6% 96.1% 



0.03 LMR 95.0% 94.7% 97.1% 96.4% 
0.03 LMA 97.0% 90.5% 98.7% 92.5% 
0.03 HM 91.1% 90.4% 96.3% 95.8% 
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Accuracy: AE bind /E bind 



Fig. 4. — Cumulative histograms for the accuracies (as defined in Eq.|6]) in our three different 
groups of models for all six metallicities in our grid. On the vertical axis, we plot the fraction 
of model data points that fall outside the accuracy quoted on the horizontal axis, hence a 
fraction of 0.1 means that 90% of our data points have that accuracy. Left panel: low- mass 
RGB, middle panel: low-mass AGB, right panel: high mass. The different line styles and 
colours indicate different metallicities, as indicated in the middle panel, and the vertical and 
horizontal dotted black lines indicate the points where the accuracy is 10% and where 90% 
of our model data points are included (10% excluded), respectively, and their intersection 
marks the point A 10 % = 90%. 



Table 4: Accuracies of our fits as in Table [3X but applied to models with stellar wind, with 
A = 1 and A = A HM (Eq.EJ). 

A = 1, A = A HM , 



Group 


Z 


A io% 


A 15% 


A io% 


A 15% 


LMR 


10~ 4 
0.02 
0.03 


98.8% 
94.3% 
93.7% 


99.3% 
96.4% 
96.4% 






LMA 


10" 4 
0.02 
0.03 


77.3% 
75.5% 
68.3% 


87.4% 
87.5% 
76.6% 






HM 


10" 4 
0.02 
0.03 


69.7% 
54.7% 
60.1% 


86.4% 
63.1% 
70.7% 


92.1% 
87.2% 
87.0% 


96.1% 
89.5% 
89.8% 
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used with arbitrary mass-loss prescriptions and still provide the accuracies listed in Table HI 
However, we find that for very strong mass loss, even this correction factor can no longer 
describe our models accurately. Our results show that when a star loses more than ~ 20% 
of its ZAMS mass, our fits become increasingly inaccurate. For each metallicity, this implies 
an upper limit in mass for which these fits can be used. We did not include the results of 
more massive models in Table H] and discuss these limitations in more detail in Sect. l4.4l 

4.4. Limits on the domain of the fits 

Here we note the limitations on the applicability of the fit functions. Outside these 
limits we do not provide fits and hence cannot estimate accuracies. 

In the absence of mass loss, Table [3] displays the accuracies of our fits for any star with 
a mass 0.8 M & < M < 100 M Q , the full range of masses we computed. However, when mass 
loss is present, we note that the parameter Ahm is only effective when stars have lost less 
than ~ 20% of their ZAMS masses. For the wind prescriptions used in this study, this limit 
can be translated to an upper mass limit of approximately 50 M Q for Z = 0.02 and Z = 0.03, 
and about 75M for Z = 10" 4 . 

In terms of radius, we have considered each model from the end of the MS to the star's 
maximum radius on the AGB, excepting only phases of radial shrinkage — most notably on 
the horizontal branch, so that any radius that may correspond to the initiation of a CE in 
a binary is covered in our study. 

The values of the metallicity considered here are discrete, and we do not provide accu- 
racies for Z < 10 -4 and Z > 0.03. Analogous to metallicities that lie between those used in 
this study, for metallicities outside this range the reader can either choose to use the value 
for Z that lies closest to the desired value, or to use an extrapolation scheme. 

5. Discussion 

In this section we discuss our results qualitatively, compare our fits to previous attempts 
to model the binding energy, note the possible effects of the choice of core definition, mixing- 
length ratio, and strong winds on the accuracy and applicability of our results. 
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5.1. Qualitative examination of the results 

Figures [2] and [3] show examples of the envelope binding energies of a selection of detailed 
stellar-evolution models (black dotted lines) compared with the results of our fits (coloured 
lines). These figures indicate that the fitting function we defined in Eq.[5] describes the 
binding energy fairly well. A full picture of the accuracy of our fits, when compared to our 
models, can be read from Tableland Fig.HJ and from Table H] and Fig.|5l for models with 
conservative evolution and models that undergo wind mass loss, respectively. 

We want to emphasise here the limitation of our results for stars that undergo substantial 
wind mass loss. When stars lose more than ~ 20% of their ZAMS mass, the binding energy 
is no longer easily related to that of a star of the same mass but without a mass-loss history. 
This is the result of an interplay between the different effects of wind mass loss on both the 
mass and structure of the envelope. Strong mass loss ultimately reduces the envelope mass, 
which decreases the binding energy compared to a conservative star of the same total mass 
and radius. However, our models indicate that the envelope structure is also affected by mass 
loss, in a way that counteracts the first effect. We find that this results in a larger binding 
energy than that of a conservative star, unless the mass loss becomes very significant. Hence, 
our compensation factor Ahm in Eq.[71 which is greater than unity, gives reasonably good 
results for small amounts of mass loss, but becomes gradually less accurate when more than 
~ 20% of the ZAMS mass of the star is lost. We note that for our choice of wind mass- loss 
prescription, this happens for stars with masses larger than > 50 M Q for Z = 0.02, 0.03 and 
for M > 75 Mq for Z = 10~ 4 , hence the most massive stars in our grids. These upper limits 
will be decreased when stronger winds are used in the stellar-evolution models. However, as 
we shall discuss in Sect. l5.3[ the dominant source of uncertainty for these high-mass stars is 
the definition of the core-envelope boundary. 



5.2. Dependence on the mixing-length ratio 



The choice of the mixing length relative to the pressure sca le height will dir e ctly in fluence 
the radius of a model star with a deep convective envelope. iDewi fc Taurisl (120001 ) looked 
at the effect of this parameter and found a weak dependence for the binding energy, with 
a higher mixing- length ratio translating into a slightly lower binding energy on the AGB. 
We compared our fits based on our standard assumption l/H p = 2.0 to a few test models 
which were computed with l/H p = 1.5 and found a similar effect: a drop in accuracy almost 
exclusive to the AGB, resulting in ~15-20% more data points outside the 10% accuracy 
range. 
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5.3. Definition of the core-envelope boundary 



The binding energy of a stellar envelope is computed by integrating the different energy 
sources (thermal, gravitational, etc.) from the core-envelope boundary to the surface of the 
star (Eq.[T]). Hence, the value of the binding energy will depend on the definition of the core 
mass M c , whi ch we chose to d e fine a s the mass coordinate at which the hydrogen abundance 
reaches 10%. iTauris fc Dewil (120011 ) consider several reasonable choices for this definition, 
one of which is identical to the definition used in our study (the central region of the star 
with X < 0.1). If we follow these authors and discard the extreme cases (first and last row) 
in their Table 1 (their Ab is almost identical to our A cnv defined in Eq.|SD, we see that while 
for low-mass stars different definitions result in similar values for M c and hence the binding 
energy, for more massive stars the results can vary quite appreciably. For a 20 M & star, the 
variation in -Ebmd can even be more than an order of magnitude. They also note that our 
choice of the definition of the core gives the lowest value for A cnv , hence the highest value for 
the binding energy, among the three options considered. 

We therefore conclude that the definition of the core mass will have little or no influence 
on our results for low-mass stars, with masses < 7 — 10 Since these s tars provide the 
majority of common-envelope donors in binaries (Ivan der Sluys et al.ll2010l ). most instances 
of CEs in nature will be little affected, as will the results for low-mass compact binaries, 
such as double white dwarfs, cataclysmic variables and low-mass X-ray binaries. On the 
other hand, for high-mass stars, the mass of the stellar core will depend strongly on how 
the core-envelope boundary is defined. In fact, in reality the exact remnant mass after the 
CE will depend on the response of both the donor st ar and the orbit to the spiral- in, an d 
hence also on the properties of the binary companion (IDeloye fc Taamll2010l ; Iivanovall2010l ). 
Therefore, for these massive stars no unique core mass can be defined and detailed stellar- 
structure and binary-evolution models are needed to determine a self-consistent outcome of 
a CE. The uncertainty thus introduced is much larger than the uncertainty we find in our 
fitting results for high-mass stars with strong winds described in Sect. 15. ll 

Finally, we would like to remark that our choice of the core-envelope-boundary definition 
is identical to that used for the fits that res ulted in the stellar-evolution prescriptions in the 
SSE and BSE codes ( Hurley et aliboool . 2002), which form the basis of a number of population- 
synthesis codes, such as StarTrack (pelczvnski et al.ll2008l ) and the latest version of SeBa 
( Nelemans et al.ll200ll ; iToonen et al.ll2010l ). Hence, even if there is no unique definition of 
the core mass, our fits will provide consistent prescriptions when implemented in the most- 
commonly used population-synthesis codes. 
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5.4. The envelope-structure parameter A e 



The purpose of this study is to provide fits of the binding energy of stellar envelopes, 
based on basic stellar parameters, so that they can be used to treat CEs in population- 
synthesis simulations, where detailed stellar-structure models are not available. So far, such 
codes have often used the so-called envelope-structure parameter, A cnv , defined by 

GMM env 



E 



bind 



(8) 



to compute the binding energy (jWebbinkl Il984t Ide Kool et al.1 119870. instead o f Eq.Q ] In 



such a case, one ne eds to assume a value for A Pnv . For example. i Nelemans et al. r|2oooh and 
Hurley et all ( 120021 ) use A cnv = 0.5, while iBelczynski et all (120081 ) choose acEA en v = 0.5, 1.0. 
Since A env has been used extensively, it is useful to give an indication here of our results 
expressed in terms of this parameter, and to compare them to the choices above. In Figure El 
we show the values of A cn v for three different stellar-evolution models and the results of 
our fits, converted to A cn v This figure alone indicates that A env is far from constant and 



will v ary f or stars of different masse s and different evolutionary phases. iDewi fc Tauris 
( 120001 ) and Ivan der Sluvs et al.l (I2006|) give more e xtensive examples of the variation of this 
parameter. In addition, Ivan der Sluys et al.l (120101 ) show that, for Z = 0.02, the assumption 
of acEAenv = 1-0 implicitly assumes that «ce > 1 for roughly 60% of the common envelopes 
in a typical stellar population, while for some cases acE > 10 is implied. These examples 
strongly indicate that a better prescription for the envelope binding energy is needed. 

The reason why we chose to develop an expression for the envelope binding energy, and 
not Acnv, may now become clear. The parameter was introduced in order to estimate the 
binding energy when it could not be computed otherwise. From now on, a reasonably accu- 
rate prescription for E^nd is available and there is no longer any need to use the parameter 



5.5. Comparison with previous work 



Previous studies that pro vided a simple prescr iption for the binding energy have been 
limited to tabulations of A env ( Dewi fc Taurisll2000l ) for stars of different masses and differ- 
ent evolu tionary stages. However, while our paper was in preparation, a paper was pub- 
lished by IXu &: Lil (120101 ) largely with the same idea as our study: to create simple and 
fast prescriptions that can provide the envelope binding energy in order to treat CEs in 
population-synthesis codes. 



The main points where our study improves on the latter paper concern the range of 
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applicability. Firstly, whereas Ku fc Lil ( 120101 ) provide fits for a selection of 14 discreet 
masses between 1M and 20 M , we use the stellar mass as a continuous parameter so 
that it can take any value between 0.8 M Q and 100 M Q . Hence, the user does not have to 
interpolate between masses, which may not be a trivial process. For example, while stars of 
1.0 M & and 1.5 M Q reach a radius of 100 R & on the RGB, a 2.0 M Q star does not. It may 
therefore be impossible to obtain the binding energy for a 1.5 M Q star at that radius on 
the RGB by interpolating the results of the 1.0 M Q and 2.0 M & models in a discrete grid. 
Secondly, they provide two choices of metallicity (Z = 0.001 and Z = 0.02), whereas we 
include six different, albeit still discreet, values, ranging from Z = 10~ 4 to Z = 0.03. We 
believe that our grid of metallicities is sufficiently dense that the nearest value can be used 
without too much loss in accuracy. Thirdly, whereas their study includes a fixed choice of 
stellar winds, our fits allow for the user's choice of wind prescription (although the accuracy 
we present is valid only within certain limits as discussed in Sect J5.1|) . 



In addition to wider a applicability, our paper presents a more extensive error analysis. 
As we have shown in Sect. 15. 1[ our accuracies are dominated by the uncertainty in wind mass 
loss for massive stars, by the uncertainty in the separation between core and envelope at the 
end of a CE, or even by the intrinsic problems in determining this separation a priori for 
stars with masses of ~ 10 — 20 M Q and up (Sect. [573]). These systematic errors are discussed 



by lXu Sz Lil (120101 ) as well. However, for the majority of common envelopes the donor mass 
is < 10 M , so that these external uncertainties are less important. For these cases the 
uncertainty in the fits will dominate, which we present explicitly in this paper. 



On the other hand, the study by Ku fc Lil (120101 ) has two advantages over ours. First, 
their work covers the HB for all masses they consider. While we include all stages of stellar 
evolution where a CE could be initiated, knowledge of the envelope binding energy on the 
HB for low-mass stars, which cannot initiate a CE, may be useful for e.g. collisions. Second, 
their prescriptions are simpler than ours, making them easier to implement. 



6. Summary and conclusions 



We provide analytical fits capable of predicting the envelope binding energy of stars 
with masses between 0.8 M & and 100 M & on the giant branches, for any given mass and 
radius, and for six discrete choices of the metallicity. These fits are based on detailed stellar- 
evolution models. In addition, we define an ad-hoc correction factor that takes into account 
wind mass loss. We provide electronic data files with the fitting c oefficients and Fortran 
routines and that are ready to use these files (jLoveridge et al.ll2010l ). 
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We find that the accuracy of our fits is better than 15% for 90% of our model data points 
in all cases and better than 10% for 90% of the data points for most cases, when fractional 
mass loss since the ZAMS is less than ~ 20%. For low-mass stars, the true accuracy will be 
close to the value we quote, whereas for high-mass stars the uncertainty in the determination 
of the core-envelope boundary will probably be the limiting factor. We conclude that our 
fits allow population-synthesis codes and other environments where detailed stellar models 
are not available to compute envelope binding energies for e.g. common envelopes quickly 
and accurately. In addition, the results presented here are more widely applicable and more 
accurate than previous prescriptions. 
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Fig. 5. — Cumulative histograms for the accuracies (as defined in Eq.[6]) in our high-mass 
models HM with wind mass loss, for a selection of metallicities in our grid. Left panel: no 
wind-correction factor applied (A = 1), right panel: with A = Ahm ; as defined in Eq.[7l 
applied. The different line styles and colours indicate different metallicities, as indicated 
in the panels, and the vertical and horizontal dotted black lines indicate the points where 
the accuracy is 10% and where 90% of our model data points are included (10% excluded), 
respectively. 



- 26 - 



t 1 1 — i — i i i i | 1 1 1 — i — i i i i | 1 1 1 — i — i i i i | 1 r 



t 

» 

t 




J I I I I ''''I I I I I ''''I I I I I ''''I I L 



1 10 100 1000 

R (R©) 

Fig. 6. — Typical values of A cnv for a selection of models in our Z = 0.02 grid. Shown 
are the LMR (solid red lines) and LMA (dashed blue lines) for low-mass models and HM 
(dash-dotted magenta lines) for the high-mass model. Masses shown are 1.04 M Q (long, high 
LMR; short, low LMA), 4.57 M (short, lower LMR; long, high LMA) and 15.3 M (HM 
track). In all cases, coloured lines indicate results from our fits, dotted, black lines show the 
detailed stellar model. 
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